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Abstract. We performed an extensive study of configuration mixing between the 
doubly excited (doorway) states and complex multiply excited states of U 28+ near its 
ionization threshold. The detailed investigation of complex spectrum and analysis of 
the statistics of eigenstate components show that the dielectronic (doorway) states 
weakly mixed with each other. However, they show substantial mixing with the 
complex multiply excited states. This situation explains the mechanism of low energy 
electron recombination with U 28+ . We calculated the energy averaged capture cross 
sections as a sum over dielectronic doorway states and found our present calculation 
interprets well the experimental recombination rates in the energy range of 1 to 100 
eV. 



PACS numbers: PACS: 34.80.Lx, 31.10,+z, 34.10,+x, 32.80.Dz 
1. Introduction 

Atomic processes are of great importance from astrophysical and other various point 
of view. Especially the low energy recombination cross sections and rate coefficients 
are required for understanding the fusion and astrophysical plasmas. In particular 
dielectronic recombination (DR) plays a significant role in astrophysics because it is 
the dominant electron-ion recombination process for most ions in low density, photo- 
ionized and electron- ionized cosmic plasmas pQ . Accurate rate coefficients are needed to 
calculate reliably the ionization balance, thermal structure and line emission of cosmic 
plasmas. Most importantly the low energy positron antiproton recombination provides 
a challenging scheme for the production of antihydrogen[2]. 

It is well-known that due to the presence of additional recombination channels 
such as DR, the rate coefficients are found to be larger than radiative recombination 
(RR) rates for many electron complex ions. In DR process the incident electron is 
captured in a doubly-excited state of the compound ion, which is then stabilized by 
photoemission. This process originally suggested by J. Sayers and was first considered 
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by Massey and Bates [3] in the study of ionospheric oxygen. Electron- ion recombination 
has been measured directly in the laboratory since early 1980's f3j. More recently the use 
of heavy-ion accelerators and electron coolers of ion storage rings has greatly advanced 
the experiment Ej. Recombination rates for various ions have been measured at 
electron energies from threshold to hundreds of electron volts (eV) with a fraction-of- 
an-eV resolution [7J |HJ UHl UH EE2 > 13J . For few-electron ions the measured rates 
were found to be in good agreement with theoretical predictions which included the 
contribution of DR resonances on top of the RR background, e.g., in He + |14j . Li-like 
C 4+ ^3] and Ar 15+ [7J^], and B-like Ar 13+ However, more complicated ions, e.g., 
Au 50+ ^1] , U 28+ PHI and Au 25+ JIB] , showed complicated resonance spectra and strongly- 
enhanced recombination rates at low electron energies. The Au 25+ ion has been studied 
in detail by Gribakin et al.|18j and Flambaum et al. [23] using statistical methods. They 
suggested that the strongly enhanced low energy electron recombination observed in 
this ion is mediated by complex multiply-excited states rather than simple dielectronic 
resonances and the dielectronic resonances play the role of doorways to the electron 
capture process. The statistical method developed by Flambaum et al. [20 a is based on 
the assumption of strong (chaotic) configuration mixing. This assumption has been 
verified by Gribakin and Sahoo JH] i n a recent study. However, for U 28+ , no theory so 
far has described the low energy DR process successfully. 

Historically the recombination rate enhancement was observed first in U 28+ |l()j. 
This measurement has been performed in a merged beam experiments at UNILAC 
accelerator in Darmstadt and at heavy-ion storage ring TSR in Heidelberg. The 
experiment found rate enhancement in the U 28+ spectrum exceeds the theoretical 
calculation by at least a factor 20 in the energy range below 10 eV. Later in 1998, they 
extended the experiment to high energies up to 420 eV and a comparison has been 
made with the theoretical calculations which is based on distorted wave approximations. 
The DR cross sections are calculated in this method are able to explain the main 
resonant features in the range 80-180 eV, but failed to identify the resonances and 
reproduce the rate at smaller energies. The cross sections involving the excitations from 
the 5s 2 5p 2 ground state configuration calculated in this method are well described by 
using either semirelativistic wave functions as found in AUTO STRUCTURE codes or 
fully relativistic wave functions in HULLAC codes. However, the resonance structure 
observed at low energies i.e. below 80 eV largely remains unexplained. They finally 
concluded that in complex ions, particularly in U 28+ ion, what are the resonances just 
above the threshold and how they contribute to the low energy recombination remain 
a 'mystery'. In this paper we performed an extensive study of the excited spectra 
and eigenstates of U 28+ near its ionization threshold and calculated recombination rate 
coefficients for electron recombination with U 28+ . We identified some of the resonances 
near the threshold those contribute significantly to the low energy recombination. The 
present results are found to be in good agreement with the experiment. This work 
develops further a statistical theory towards the full understanding of the mechanism 
of low energy electron recombination with U 28+ and other similar complex ions. 
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2. Many-electron excitations 

We consider the recombination of an electron with U 28+ . Due to electron correlation the 
slow electron can be captured in one of the excited states of the compound ion U 27+ . 

It has 65 electrons and its ground state configuration belongs to Is 2 5s 2 5p 3 . Figure 

1 shows the spectrum of relativistic nlj orbitals obtained by relativistic Hartree-Fock 
calculations. Atomic units (a.u.) are used unless otherwise stated. All the occupied 
orbitals below fermi level (~ 31.19 a.u.) are obtained by a self consistent calculation 
of U 28+ ground state. Each of excited orbitals above the fermi level are calculated by 

placing one electron on to it in the frozen Is 2 5s 2 5p 2 (U 28+ ) core. Our configuration 

interaction (CI) calculation shows that the ground state of U 28+ and U 27+ ions are 
characterized by their total angular momentum J = 0.0 and 1.5 respectively. The 
difference between their total energies —27741.40 a.u. and —27771.37 a.u., gives us an 
estimation of ionization threshold (I.T.) = 29.97 a.u. (815.184 eV). 

The excited states of the ion are generated by transferring one, two, three, etc., 
electrons from the ground state into the empty orbitals above fermi level. Since we 
are interested in the excited spectrum of the ion near its ionization threshold (29.97 
a.u.), we consider it as a system of having only 29 electrons above the frozen Kr-like 
ls 2 ...Ap 6 core. The number of many-electron states obtained by distributing 29 electrons 
over 40 relativistic orbitals, 4g? 3 / 2 through to 8(79/2 are huge in number. It is practically 
impossible to perform a CI calculation for all of them. We construct the excite spectrum 
by using mean field approach by calculating their mean energies E iy and number of many 
electron states Aj associated with each of them: 

Ei = E core + e a n a + ^ — TTT, ab > ^ 

a 

gJ 



1 + $ab 

a<b 



nj(9a ~n a )V 



(2) 



where n a are the orbital occupation numbers of the relativistic orbitals in a given 
configuration and J2 a n a = n - e a =< d\H core \a > is the s 

ingle-particle energy of the orbital a in the field of the core, g a = 2j a + 1, and U a b 
are the average Coulomb matrix elements for the electrons in orbitals a and b ( direct 
minus exchange): 



U, 



ah 



d(0) _ V^r n(A) \ja }b A 
n-abab / j °P IX abba \ 1 _1 n 
A L 2 2 



(3) 



9a ~ Oab 

R-abba * s ^ ne two-body radial Coulomb integral of A multipole, and 5 P = 1 when 
l a + h + A is even and otherwise. Using Eqs. (l)-(3) we obtained about 353 
configurations within ± 1 a.u. of ionization threshold. They comprise a total of 
1.9 x 10 5 many electron states. The single-particle spectrum of U 28+ does not show large 
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gaps. Owing to the "gapless" single-particle spectrum, the density increases rapidly as 
a function of energy, as described by the Fermi-gas-model ansatz [2~T] 

p{E) = AE~ u exp{a^), (4) 

with A = 0.0885, v = 2.33, and a = 3.99 a.u. [23], where E is the energy above the 
ground state in atomic units. Figure 2 shows the level density calculated by averaging 
with a Gaussian with 1 a.u. variance. This figure depicts the level densities for both odd 
and even parity configurations. Also included in the figure are only the doubly excited 
configurations of even and odd parity found within ± 1 a.u. of ionization threshold. It 
is found that the level density of the odd configurations near its ionization threshold is 
of the order is of 1 x 10 5 and 
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Figure 1. Level densities in U 28+ within ±1 a.u. of Ionization threshold. Thick solid 
line: odd configurations, Thin solid line: odd dielectronic configurations, Dash-dot 
line: even configurations and dashed line: even dieclectronic configurations. 

those of even configurations it is of the order of 1.3 x 10 4 . This provides an evidence 
that this system is characterized by huge level density. This system may be compared 
with Au 25+ , in which the level density near its ionization threshold is about 1 x 10 7 |18j. 
This shows that U 28+ ion is less complicated than Au 25+ from the view point of dense 
excited spectrum. Apart from this one can estimate the mean level spacing from the 
distribution of total angular momentum J as shown in Figure 3. 

Ionic eigenstates are characterized by their total angular momentum and parity J*", 
and are 2 J + 1 times degenerate. Therefore the total level density can be broken into 
a sum of partial level densities: p(E) = ^j^(2J+ l)pj^(E). The excitation spectrum 
of U 28+ near the ionization threshold, E = I « 29.97 a.u., contains many J ranging 
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from i to y. Their distribution is in agreement with statistical theory fU\ I22|. which 
predicts that at a given energy pj* are proportional to the function 
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where J m is the most abundant J value. Numerically for U 27+ we find J, 
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Figure 2. Densities of states with different J near ionization threshold of U 28+ . 
Histograms show their distribution and the solid line is non linear fitting obtained 
from Eq.(5) which predicts the most abundant value of J ~ 3.2 and the proportionality 
constant ± 11075.9. 

Using Eq. (jSJ) one can estimate the partial densities by pj = pj+ + pj- ~ 
f(J)p/ (2J+ 1), where (2 J + 1) is an average over f(J). For the most abundant angular 
momenta J ~ J m , and assuming pj+ m p 3 -, we have pj^(E) = AjnE~ u exp(av^E'). Near 
the ionization threshold this gives pj* ~ 1.8 x 10 3 a.u. (Fig. 3), which means that the 
spacing between the multiply-excited states with a given J n is small: D = l/pj* ~ 15 
meV. Whereas in Au 25+ it is about 1 meV, i.e. 15 times larger than that found in Au 25+ . 
This situation explains why individual resonances appear in recombination rates of U 28+ 
but the same is not observed in Au 25+ experimentally even at an energy resolution of 0.1 
eV. However, the large density of multiply-excited states is only a "kinematic" reason 
behind the experimental observation. To explain the fact that the electron can actually 
be captured into these states, we need to analyze the dynamics of electron capture and 
show that the residual Coulomb interaction between the electrons (i.e. that beyond 
the mean field) makes for an efficient capture and accounts for the observed enhanced 
recombination rate. 
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3. Configuration mixing 

Taking into account the fact that the rasidual Coulomb interaction is the key problem 
in many-electron processes, we construct the basis of many-electron states from 
single-particle (e.g., Dirac-Fock) orbitals, and solve the eigenvalue problem for the 
Hamiltonian matrix = (<f>i\H\<f>k), which yields the eigenvalues E v and eigenstates 
\^ v ) = Ylk^k l^fc) °f ^ ne system (configuration interaction method). We performed 
two sets of model calculations. One includs all the configurations within ±1 a.u. of 
ionization threshold of U 28+ , which produces 2516 states for J^—l sequence. Similarly 
the second calculation includes all the dielectronic (doorway) states within ±1 a.u. of 
ionization threshold and produces 108 states for J 7r =| sequence. This shows that the 
number of states associated with the dielectronic configurations for a given J n are not 
large in number. As a result we performed the full CI calculations for Hamiltonian 
matrix of size 108 and 2516 respectively and obtained the eigen values and eigenstate 
components. To study the mixing between the doubly excited states we analyze the 
eigenstate components by calculating the weight of a given doorway configuration shown 
in Figure 4(a)-(e). The weight (w) of a doorway state can be calculated as £fc=i,jv c |Cfc| 2 j 
where N c is the number of states in each configuration. It has been shown in Ref 19 
that when there is strong and uniform configuration mixing the weights significantly 
reduce from 1, but in the present case the weights go down from 1 but not significantly. 
So one can say that these doubly excited (doorways) states weakly mix with each other. 
When they are included in the large calculation (2516 states), their weights significantly 
go down from 1 as shown in Figure 4(f)-(j). This gives a signature of strong mixing i.e., 
the doorways show a significant mixing with multiply excited states. It has also been 
observed that the multiply excited configurations mix with each other quite comfortably. 
This mixing is mainly responsible for an enhancement of recombination rates over RR. 
In the recombination cross sections, the doorways which do not mix completely with 
either multiply excited states or with other doorways (in a sense they remain 'isolated') 
appear in the form of narrow peaks (resonances) which explains the experiment as well 
as the theory [T7j. 

On the other hand, when the level density is high and the two-body interaction 
is sufficiently strong the system is driven into a regime of many-body quantum chaos, 
where the effect of configuration mixing can be described statistically in the case of 
Au 25+ ^Hj- This regime is characterized by the following |2"o] 127]. (i) Each eigenstate 
contains a large number iV of principal components ~ 1/ \/N, corresponding to the 
basis states which are strongly mixed together, (ii) Owing to the strong mixing, the 
only good quantum numbers that can be used to classify the eigenstates, are the exactly 
conserved total angular momentum and parity J n and the energy, (iii) The degree of 
mixing in this regime is in some sense complete, i.e. all basis states that can be mixed 
(within a certain energy rage, see below) are mixed. The notion of configurations based 
on the single-particle orbitals becomes largely irrelevant for the purpose of classifying 
the eigenstates. Each eigenstate contains substantial contributions of a few nearby 
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configurations. As mentioned above, in the present case we found that there is a weak 
configuration mixing between the dielectronic doorway configurations and complicated 
multiply excited states, and the multiply excited states show a substantially strong 
configuration mixing with each other. But the degree of mixing is not sufficiently strong 
to drive the system into chaotic regime. 
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Figure 3. Weights of odd dielectronic doorway states obtained from two different 
calculations. Figure (a)-(e) are the weights of doorways obtained from a calculation 
involving 108 states. Figure (f)-(j) are the weights of the same doorways obtained from 
a calculation involving 2516 states 

(a) :4^ /2 4^ /2 4/| /2 4/ 7 7 /2 5^ /2 5p2 /2 5pi /2 7/ 5 1 /2 (7V c =4) 

(b) :44 /2 4^ /2 4/| /2 4/ 7 7 /2 5^ /2 5p2 /2 5pi /2 7/ 7 1 /2 (iV c =4) 

(c) :44 /2 4^ /2 4/ 5 6 /2 4/ 7 8 /2 5^ /2 5pi /2 5^ /2 6pi /2 (iV c =2) 

(d) :4^ /2 4^ /2 4/ 5 6 /2 4/ 7 8 /2 5 sr/2 5^ /2 54 /2 7/ 5 1 /2 (iV c =2) 

(e) :44 /2 4^ /2 4^ /2 4/ 7 8 /2 5^ /2 5^ /2 5 /5 i /2 6p} /2 (iV c =2) 

N c is the number of states associated with each of the configurations. 

This feature can be studied from the inverse participation ratio (IPR): 
£j=i,jVc|Cfc |~ 4 - Figure 5 shows the inverse participation ratio for the two sets of model 
calculations. The top figure depicts the inverse participation ratio of the doorways 
only and the middle one for the large (2516 states) calculation that includes all the 
configurations (both dielectronic and multiply excited states). We found most of the 
doorways are associated with number of states either 4 or 2 or even less. The top figure 
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Figure 4. Inverse participation ratio (IPR). Top figure: for doorways having 108 
states, the middle figure: for all configurations having 2516 states and the bottom one 
is same as middle figure but in an extended scale. 



shows that the inverse participation ratio lies flat between 2 and 4 except at the energy 
range where the eigen energies are close to the ionization threshold. This indicates 
that a few number of the doorways having energies close to the threshold participate in 
mixing. However the middle figure shows a picture of strong but non uniform mixing 
which involves a lots of multiply excited states. When we included the IPR of doorways 
with the IPR obtained from large calculation as shown in the bottom figure, it lies well 
below as indicated by solid circles. It also provides us with an information that some 
of the doorways take part in mixing with multiply excited configurations and other 
doorways are weakly mixed or even remain isolated. 

This strong mixing takes place in a certain energy range \E^ — E v \ < r spr , where 
Ek = Hkk is the mean energy of the basis state and r spr is the so-called spreading width. 
More precisely, the mean-squared value of as a function of E k — E u , can be described 
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by a Breit-Wigner (BW) formula 
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Figure 5. Spreading of doorway configurations in complicated multiply excited states 
obtained from 251 6 x 2516 calculations for J 71 " = 3.5~ sequence. Each of the doorway 
configuration spreads as (a) 0.1433 a.u., (b) 0.1310 a.u., (c) 0.1661 a.u., (d) 0.2321 a.u. 
and (e) 0.2310 a.u.. These doorways are same as indicated in Fig4. 



The mean-squared components are obtained by averaging over the basis states 
associated with each of the doorway configurations and are plotted as a function of 
eigen energies. We obtained r spr of doorways (close to the ionization threshold) from 
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the BW fit which is shown in Figure 6(a)-(e). From the BW fit we observed that the 
r spr is not constant and it varies from 0.1 a.u. to 0.2 a.u. in these doorways. There may 
be doorways which even show much less spreading. Roughly one can say the mixing 
takes place within 0.1 a.u.. So we calculated the recombination rates with two different 
values of spreading width which can be found in section IV. It may be pointed out that 
though the single BW fit is not as accurate as it should be, still it gives an estimation 
of the important quantity i.e. r spr . It is worth mentioning that the spreading width 
found in Au 25+ is about 0.5 a.u.^HJ and these value does not change if one performs 
calculations by adding more configurations [T§] . 



4. Recombination 



4-1. Theory 

For low-energy electrons the contribution of the autoionising states (resonances) to the 
recombination cross section is given by (see, e.g., Ref. [24] ) 

_ n v 2J+1 T ( M a) 
° r fe2^2(2Ji + l) (£-<) 2 + r2/4' U 

where e = k 2 /2 is the electron energy, Jj is the angular momentum of the initial (ground) 
target state, J are the angular momenta of the resonances, e v — E v — I is the position of 
the z/th resonance relative to the ionization threshold of the compound (final-state) ion, 
and ri , and T u = T„ + V„ are its autoionisation, radiative, and total widths, 
respectively j2Sj- When the resonance spectrum is dense, a r can be averaged over an 
energy interval Ae which contains many resonances, D <C Ae <C e, yielding 

_ _2vr 2 ^ 2J + 1 / rj r T| a) \ 

° T k 2 2^ 2(2j i + i)z> \rw + it7 ' [ ) 

where (. . .) means averaging. If the fluorescence yield, Uf = / (Tv +T^), fluctuates 
weakly from resonance to resonance (see below), one can write a r = a c ujf, where 



7T 



2 



a. 



y- ( 2J + 1 ) r(a) (9) 



k 2 ^ (2J t + l)D 



is the energy-averaged capture cross section, and is the average autoionisation width. 

In a situation when there is a strong configuration mixing between the dielectronic 
doorway states and multiply excited states, the capture cross sections can be obtained 
as a sum over the single-electron excited states a, j3 and hole states 7, as well as the 
partial waves Ij of the continuous-spectrum electron e. As a result, we have 

tt 2 \- T spr x - (a,/3\\V x \\j,elj) 



Or 



(a,/3||y A || 7 ,di) - (2A+ l)J2(-l) X+X ' +1 { i 3a 3 } (a,f3\\V x ,\\elj,d) 
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where e a , £p and e 7 are the orbital energies, the two terms in square brackets represent 
the direct and exchange contributions, and (a, /3\\V\\\y, elj) is the reduced Coulomb 
matrix element (see Ref.|20j). 

It is assumed that the energies of dielectronic doorway states relative to the 
threshold is given by e a + Ep - e t A more accurate value can be obtained by using mean 
field energies (configuration energies) of doorway configurations in Eq.(lO). The effect of 
using these two different energies can be found in next subsection. We have also shown 
the the qualitative difference between the results obtained by using two different values 
of T spr . Because r spr is well defined in the case of a strong and chaotic configuration 
mixing. However, it can not be properly defined if the mixing of the configurations is 
weak and non-uniform. 

Equation (fTU|) is directly applicable to targets with closed-shell ground states. If 
the target ground state contains partially occupied orbitals, a factor 



where n a , np, and n 7 are the orbital occupation numbers in the ground state $«, must 
be introduced on the right-hand side of Eq. (fTUj). Steps similar to those that lead to Eq. 
()1L)|) were used to obtain mean-squared matrix elements of operators between chaotic 
many-body states [23 12H] ■ 

The chaotic nature of the multiply-excited states can also be employed to 

(r) 

estimate their radiative widths T v ■ Electron-photon interaction is described by a single- 
particle dipole operator d. Any excited electron in ^> v may emit a photon, thus leading to 
radiative stabilization of this state. The total photo-emission rate Tu can be estimated 
as a weighted sum of the single-particle rates, 



and /3, and (. . .)„ is the mean occupation number factor. Since ^> u have large numbers of 
principal components N, their radiative widths display small 1/y/N fluctuations. This 
can also be seen if one recalls that a chaotic multiply-excited state is coupled by photo- 
emission to many lower-lying states, and the total radiative width is the sum of a large 
number of (strongly fluctuating) partial widths. A similar effect is known in compound 
nucleus resonances in low-energy neutron scattering [21J. 

There is a certain similarity between Eqs. fllOl) and ()12)) and those for autoionisation 
and radiative rates obtained in a so-called configuration-average approximation [29\. 
In both cases the answers involve squares or products of two-body Coulomb matrix 
elements [see the direct and exchange terms in Eq. (jlO))] . or single-particle dipole 
amplitudes [Eq. ([12p]. However, there are a number of important differences between 
the present results and the configuration-average approximation. The latter considers 
dielectronic recombination and introduces averaging over configurations as a means of 
simplifying the calculation. The DR cross section is averaged over an arbitrary energy 




(11) 




where u)p a = Ep 



— E a > 0, (a;||d||/3) is the reduced dipole operator between the orbitals a 
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interval Ae, and only the configurations within this energy range contribute to the 
average. Effects of configuration mixing as well as level mixing within a configuration 
are neglected. 

It is important to compare the radiative and autoionisation widths of chaotic 
multiply-excited states. Equation (fT2j) shows that is comparable to the single- 
particle radiative widths. On the other hand, the autoionisation width T^ a \ is suppressed 

M 2 

by a factor ~ iV -1 relative to that of a typical dielectronic resonance. Therefore, 

in systems with dense spectra of chaotic multiply-excited states the autoionisation 
widths are small. Physically this happens because the coupling strength of a two- 
electron doorway state to the continuum is shared between many complicated multiply- 
excited eigenstates. As a result, the radiative width may dominate the total width of 
the resonances, ^> T^ a \ making their fluorescence yield close to unity. However, 
in the case of U 28+ , T^' ~ r ■ Our numerical results for the recombination of U 28+ 
presented in Sec. I4.2[ confirm this scenario. 

The resonance recombination cross section should be compared with the direct 
radiative recombination cross section 

32tt Zf. ( Zi\ , 1Q , 
^ = -^Trln -T , (13) 



3V3c 3 k 2 \n k 

obtained from the Kramers formula by summing over the principal quantum number 
of the final state JH]- Here Zj is the ionic charge {Z^ = 28 for U 28+ ), and Uq is the 
principal quantum number of the lowest unoccupied ionic orbital (n = 5). Note that 
the direct and energy-averaged resonance recombination cross sections of Eqs. ()13|) and 
(JHJ) have similar energy dependences. 



4-2. Numerical results 

Numerical calculations of the cross section from Eqs. (|10|) and involve summation 
over the orbitals shown in Fig. 1 with electron partial waves up to in/2- The 
results of the calculations for the recombination rates are displayed in Figure 7. We 
calculated the recombination rates using two different energies, i. e., configuration 
energies and relativistic energies. It is found that both the calculations are in good 
agreement. However, the magnitude of recombination rates obtained from both the 
present calculations lie above the experimental rates. This is due to the fact that we used 
ujf — 1 in the present calculations similar to that in Au 25+ . It is clear that in the present 
case Uf is less than 1 and remains constant throughout the energy range considered. 
The exchange contributions are found to be about 200 times smaller than the direct 
one so we did not include these when calculated the final rates. The radiative rates are 
found to be smaller in magnitude in comparison to DR rates and have almost negligible 
interference with the resonances. It may be pointed out that the present calculation with 
spreading width T spr ~ 0.15 a.u, does not show enough resonance structures as have been 
observed in the experiment. Hence, we performed another calculation with T spr ~ 0.05 
a.u. and compared with the results obtained using spreading with 0.15 a.u. shown in 
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Figure 8. It is clear that the results do not show any change in magnitude. However, 
the calculation with T spr ~ 0.05 interprets the resonance peaks very well which are in 
reasonable agreement with the experiment, though the position of peaks are different 
from the experiment. This is because the present energies of the doorway states are 
approximate and a few eV relative error is expected in this approach. The presence of 
narrow peaks in rate coefficients can be interpreted as: the dielectronic states which 
play the role doorway to the electron capture process weakly mix with each other as 
has been discussed in previous section and appear as single peaks. Because in statistical 
calculations one peak corresponds to one doorway state. It may be recalled that the 
distorted wave calculation jTJj, predicts the recombination rates in agreement with the 
experiment above 80 eV and below this energy the results are smaller in magnitude. 




10 20 30 40 50 60 70 80 90 100 
Electron energy (eV) 



Figure 6. Electron recombination rate on U 28+ . Open circles connected by solid line: 
Calculation 2 (using configuration energy and T sprea( i = 0.15 a.u.), Dot-dashed line: 
Calculation 1 (using HFD energy and T spreae i = 0.15 a.u.), Solid lines : Experimental 
data, Solid line connected by plus sign : Exchange contribution (using HFD energy 
and T spreac i = 0.15 a.u.), Dotted line: Radiative rate. 

Table ^ lists the most important dielectronic doorway contributions to the 
dimensionless sum in Eq. f!10|) . which also determines the ratio of the autoionisation 
width to the spacing between the resonances. In total they account for about two thirds 
of the total cross section. Although these transitions have been selected according to 
the size of their contribution, their energies from two different calculations are close 
to the threshold, as seen in the last two columns in Table H Indeed, the spreading of 
configurations discussed in Sec. |2] allows configurations near the threshold, \AE\ < r spr , 
to contribute. On the other hand, the contribution of configurations lying far away 
from threshold, \AE\ ^> r spr , is suppressed. It may be pointed out that the theory 
of Mitnik et. al. considers the excitations from 5s orbital only. However, the 
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Table 1. Electron orbitals which give the leading contribution to the low-energy 
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[l]a and (3 are the 



excited electron orbitals, and 7 is the ground-state hole of the dielectronic doorway 
state; elj is the partial wave of the incident electron. [2]Direct term contributions to 
the dimensionless sum in Eq. (110)1 . a c k 2 /7t 2 , using spreading width = 0.15 a.u. and 
configuration energy , with magnitudes greater than 5 x 10~ 3 . [3]A£? = e a + Bp — e 7 
is the mean-field energy of the doorway state relative to the threshold. [4]A£J is the 
configuration energy relative to the threshold. 
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present calculations show that excitations involving 4/ orbitals are more important 
especially at threshold since they contribute significantly to the recombination cross 
sections. As expected, the energy dependences of the resonance and direct recombination 
rates are very close, although the latter is about 200 times smaller. There is a good 
overall agreement between the resonance rate and experimental data at electron energies 
between 1 eV and 100 eV. To compare with Au 25+ , we calculated the rate coefficients for 




50 100 
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Figure 7. Electron recombination rate on U 28+ . Open circles connected by solid 
line: Calculation 2 (using configuration energy and T spreae i — .15 a.u.), Solid line: 
Calculation 2 (using configuration energy and r spreo d = 0.05 a.u.), Dense solid line : 
Experiment, and Dotted line: Radiative rate. 

the electron energy ranging from 1 eV to 100 eV. Because the previous calculation 20 
was restricted to the energy range below 1 eV. Figure 9 shows results three different 
calculations and the experimental data. It may be seen that the present (we call our 
best calculation) calculation uses configuration energies shows a very good qualitative 
agreement with the experiment in comparison to the other two calculations. The RR 
is found be of smaller in magnitude. The most striking feature in this graph is that 
the experimental data quickly departs from the theoretical values as one goes to higher 
energies. This can be interpreted as that in this energy range the fluorescence yield 
certainly goes down from 1 which can be well understood from the Figure 10. In this 
figure we show the configuration energies as a function of total angular momentum J. 
It shows that due to the presence of a lots levels, many inelastic channels open up and 
they participate strongly. This situation makes the fluorescence yield less than 1 and 
hence the theory predicts recombination rates higher in magnitude. 

Which in turn just opposite to the case of U 28+ as shown in Figure 11. 

Finally, we compare in table II, the present energies of the lowest configurations 
in U 28+ with the energies calculated by using HULLAC code^Zj. As mentioned 
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Figure 8. Electron recombination rate on Au 25+ . Open circles connected by solid 
line: Calculation 2 (using configuration energy and T sprea d = 0.5 a.u., Thick Solid line 
(using HFD energy and T S p reac i — 0.5 a.u.), Thick dotted line: Original calculation (as 
in PRA 66, 012713 (2002)), Dense solid line : Experiment, and Dotted line : Radiative 
rate. 
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Figure 9. Spectrum of low lying excited states of Au 25+ . All the levels shown belong 
to the 41* configurations. The ground state is characterized by J=6. 



earlier that in RefjTTj, energies have been calculated by using different codes such 
as multiconfiguration Hartree-Fock (MCHF), the AUTOSTRUCTURE code in the 
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Figure 10. Spectrum of low lying excited states of U 28+ . All the levels shown 
belong to 5s 2 5p 2 , 5s5p 3 and 5s 2 5p5d configurations (see table II). Its ground state 
is characterized by J=0. 



perturbative-relativistic [AS (PR)] and in the semirelativistic [AS(SR)] mode and the 
HULLAC code. On comparison we found the present energies are close to those 
obtained by Mitnik et. al. J7j using HULLAC code since both the calculations are 
fully relativistic. 



5. Summary and outlook 

We have shown that the dielectronic states weakly mix with each other and show 
a substantial mixing with complicated multiply excited states. This explains the 
mechanism of low energy recombination of electron with U 28+ . The present results are 
found to be in good agreement with the experimental data. We found that one must not 
ignore the excitations from 4/ orbitals while considering the low energy recombination 
of U 28+ ion. On detailed study, we found the configuration mixing between the doubly 
excited states and multiply excited states is not complete or uniform. Although this 
theory is valid in a system having an extreme degree of configuration mixing |20j, 
still it predicts quite good results for other systems where there is not enough strong 
configuration mixing, for example U 28+ . This work develops further a statistical theory 
to study the low energy recombination process. 
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Table 2. Energies of the lowest configurations in U 28+ . The present calculated energies 
obtained by using Hartree Fock Dirac (HFD) code are compared with the energies 
calculated using HULLAC code [18]. The energies are in eV. 
Configurations Term(J) HULLAC HFD (present) 
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